lambda = -6000000
list = 0:0.001:3;
plot(list, exp(lambda*list)*(1.5 - 1) + cos(list), '-r');
print(gcf, '-depsc','IVP_plot.eps')